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ABSTRACT 

Radio point sources are strongly correlated with clusters of galaxies, so a significant fraction of the thermal 
Sunyaev-Zel'dovich (SZ) effect signal could be affected by point source contamination. Based on empirical 
estimates of the radio galaxy population, it is shown that the rms temperature fluctuations of the thermal SZ effect 
could be underestimated by as much as 30% at an observing frequency of 30 GHz at £ > 1000. The effect is larger 
at higher multipoles. If the recent report of excess power at small angular scales is to be explained by the thermal 
SZ effect, then radio point sources at an observing frequency of 30 GHz must be a surprisingly weak contaminant 
of the SZ effect for low-mass clusters. 



1. INTRODUCTION 

Recent and upcoming measurements of anisotropics in 
the cosmic microwave background (CMB) may be sensitive 
enough (a few fiK) and at high enough angular resolution (a 
few arcminutes) to detect fluctuations in the microwave back- 
ground due to unresolved distant and faint clusters of galax- 
ies, due to the thermal Sunyaev-Zel'dovich (SZ) effect. As the 
largest expected signal of "secondary" anisotropies (anisotropy 
arising from processes at low redshift), such a detection will 
mark an important milestone, as well as provide important con- 
straints on cosmological parameters and the thermal history of 
intra-cluster gas. 

Significant work has gone into the problem of extracting 
the signal of the "primary" anisotropies (those imprinted at 
Z ~ 1100) from the expected foreground sources of contam- 
ination, such as dust, star-forming galaxies, and radio point 
sources. A major aid in this extraction is the lack of correla- 
tion between the signal and the contaminant. In searching for 
the signal of secondary anisotropies, this is no longer the case. 
While this does not provide any deep or profound problems, 
it does provide reason for caution in interpretation of detected 
signals in the presence of foregrounds. 

It is well-known that clusters of galaxies observed at frequen- 
cies near 30 GHz (wavelength 1 cm) are likely to contain ra- 
dio point sources (Birkinshaw 1999), primarily due to radio 
emission from the cluster galaxies themselves (Cooray et al. 
1998). Therefore, the signal of the thermal SZ effect can be 
significantly diluted, depending on how these point sources are 
handled. There are two possible modes of dilution: removal 
of point sources can lead to removal of SZ signal, and point 
sources that are not removed can fill in an SZ decrement. In 
this work we estimate the expected dilution of the thermal SZ 
effect signal at an observing frequency of 30 GHz, appropriate 
for an experiment such as CBI (Padin et al. 2000). 

While the focus in this paper will be contamination of the 
thermal SZ effect by radio point sources, at higher frequencies 
dusty starbursting galaxies could be important contaminants of 
the SZ signal (Blain 1998), due to gravitational lensing ef- 
fects. These same contaminants will be troublesome for sig- 
nals of lensing in the CMB. Multifrequency observations al- 
leviate these concerns somewhat, but the uncertain and het- 
erogeneous spectral behavior of both radio and submm point 
sources present a challenge for precise measurements of sec- 



ondary anisotropies. 

The thermal SZ effect and its angular power spectrum is re- 
viewed in §2. In §3 we outline the relevant statistics of radio 
point sources in galaxy clusters and apply them to the thermal 
SZ power spectrum in §4. The implications of these results are 
discussed in §5. 

2. THERMAL SZ FROM GALAXY CLUSTERS 

The thermal SZ effect (Sunyaev and Zel'dovich 1972; 
Birkinshaw 1999; Carlstrom et al. 2000) arises from Compton 
scattering of cool CMB photons with hot electrons in the deep 
potential wells of galaxy clusters. A striking feature of the SZ 
effect is its unique spectrum. Relative to the CMB, the SZ ef- 
fect manifests itself as a deficit of photons at frequencies below 
about 218 GHz and as an excess at higher frequencies. The 
decrement at low frequencies is particularly useful, since there 
are very few astronomical signals which show up as "holes" in 
the sky. 

The temperature decrement (or increment), ignoring rela- 
tivistic corrections (Rephaeli 1995), is given by 

^^■ = f(x)y = f{x)[n e ^^a T d£, (1) 
Tcmb J m e c z 

where y is the Compton y parameter, n e is the electron num- 
ber density, T e is the electron temperature, and x = is the 
observing frequency in natural units. The frequency dependent 
factor is given by f(x) = (*frz[-4). When integrated over the 
entire angular extent of the cluster, the SZ flux has a scaling 
of S sz oc d~^ 2 f g hMT e , with the angular diameter distance in 
units of /r'Mpc, f g the gas mass fraction, h the Hubble constant 
in units of 100 kms" I Mpc" 1 and M is the cluster mass in units 
of h~ l MQ. Assuming the virial relation T oc M 2 / 3 gives the SZ 
flux scaling with mass as M 5 / 3 . 

The primary anisotropies in the CMB are thought to be 
Gaussian in nature, and can therefore be characterized entirely 
through the angular power spectrum. This is not true for ei- 
ther foregrounds or secondary anisotropies. The angular power 
spectrum is therefore not nearly as useful, but it does present 
a common language with which to work and we will adopt it 
here. 

Adopting the flat-sky approximation, the usual multipole ex- 
pansion becomes a two dimensional Fourier transform, with 
£ = 2%R UV , where R uv = Vu 2 + v 2 is the radial distance in the 
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Fourier plane and u and v are the Fourier conjugate variables 
to 9 X (e.g., R.A.) and 6 y (e.g., Dec.) on the sky. The variance 
of the Fourier amplitudes at radius R m , is equal to q. Denoting 
the Fourier transform of the cluster SZ temperature decrement 
profile as T, the angular power spectrum of a galaxy cluster 
at the center of a field is simply q = T 2 , assuming azimuthal 
symmetry for the cluster. 

Spatial correlations between galaxy clusters are negligible 
for t > 100 (Komatsu and Kitayama 1999) and in this work we 
are interested in I > 1000, so the correlations can be safely ne- 
glected. The angular power spectrum can therefore be thought 
of as Poisson shot noise, where each "shot" has an angular pro- 
file. For a Poisson process, at each position in the Fourier plane 
a randomly placed source will have an amplitude T but will 
have a random phase. The collection of sources will thus con- 
stitute a random walk of the Fourier amplitude with step size 
T, leading to a final q = NT 2 , where N is the number of galaxy 
clusters. 

As an integral over redshift and cluster mass, this can be ex- 
pressed as 

(2) 

where cIa(z) is the angular diameter distance, H(z) is the Hub- 
ble constant, and dn/d InM is the differential comoving number 
density per log interval in mass. A very fast method to obtain f 
for somewhat realistic cluster profiles with azimuthal symme- 
try is to take advantage of fast Hankel transform routines that 
exist (Anderson 1982). Note that here, and everywhere below, 
we work in units where q has units of [iK 2 . 

The power spectrum calculation exactly follows the proce- 
dure of Holder and Carlstrom (2001) . We choose cosmologi- 
cal parameters Q m = 0.3, f^A = 0.7, a% = \,h = 0.7, flbh 2 = 0.02, 
n = 1, and zero neutrino mass. The differential comoving num- 
ber density is adopted from recent fits to large numerical sim- 
ulations of structure formation (Jenkins et al. 2001), and is 
a function of the variance on mass scale M. This variance was 
calculated using the power spectrum for our adopted cosmolog- 
ical model derived from the fitting functions of Eisenstein and 
Hu (1999) . 

As a cluster model, we adopt the simple toy model of Holder 
and Carlstrom (2001), with a density profile of the form n e oc 
l/(r 2 +R 2 ), with r c a core radius and R v the virial radius, de- 
rived from the spherical collapse model (Lahav et al. 1991). 
The relation between core radius and virial radius is taken to 
be a constant value (ten), roughly as would be expected for the 
case of self-similar evolution of the cluster population. The 
gas temperature as a function of mass was taken from hydro- 
dynamical simulations (Bryan and Norman 1998) and the gas 
was assumed to be isothermal. 

For this work, we are interested in the relative effects of ra- 
dio point sources on the SZ angular power spectrum, so the 
details of the power spectrum are not crucial. We have ver- 
ified that the results below are robust to the choice of model 
for generating the SZ power spectrum. The angular power 
spectrum resulting from our recipe is in broad agreement with 
results from large cosmological hydrodynamical simulations 
(Springel, White, and Hernquist 2000). Temperature gradients 
will affect the details of the shape of the peak of the power spec- 
trum, while shifting the normalization of the mass temperature 
relation introduces a direct scaling of the amplitude of the angu- 
lar power spectrum. For example, using the observed normal- 



ization of the mass-temperature relation (Finoguenov, Reiprich, 
and Bohringer 2001) leads to an rms temperature fluctuation 
that is approximately 50% higher. Given the uncertain relation 
between the observed X-ray temperatures and the (most rele- 
vant for our purposes) mean electron temperature, we choose 
to use the normalization from simulations. 

Putting the pieces together, the integrand of Equation 2 is 
shown as a function of mass and redshift in Figure 1. Self- 
similar evolution of the intra-cluster medium was assumed, but 
the gas evolution history has very little effect on such a plot. At 
I = 1000, most of the signal is coming from somewhat massive, 
relatively nearby clusters, while the dominant contribution to 
higher I is coming from distant low-mass clusters, with a sig- 
nificant tail extending to z = 2. 

It is easy to see why radio point sources could be a prob- 
lem. At I = 1000, the signal is coming from nearby clusters, 
where each point source should be relatively bright, and from 
fairly massive clusters, which have more cluster members and 
therefore could be expected to have more point sources. These 
clusters also have stronger SZ emission and are more extended 
on the sky, which will somewhat mitigate the point source con- 
tamination. At higher £, the clusters have fewer point sources 
and the ones that they do have are diluted by the luminosity 
distance. However, the cluster signal is weaker, making point 
sources relatively more important. Thus, over the whole range 
of the peak of the SZ angular spectrum it is expected that radio 
point source contamination might be important for experiments 
at low frequencies. 
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FIG. 1. — Relative differential contribution to the angular power spectrum, 
dct/dlnMdz, normalized to peak contribution. Solid contours show contribu- 
tions to £ = 1000, with contour levels of 10% of the peak contribution. Dotted 
contours show contributions to £ = 4000, also with contour levels of 10% of 
peak. 

3. RADIO POINT SOURCES IN GALAXY CLUSTERS 

Radio point sources are well measured and catalogued in 
terms of their flux distribution and source density at an ob- 
serving frequency of 1.4 GHz, due to the NVSS (Condon et al. 
1998) and FIRST (White et al. 1997) surveys. Unfortunately, 
radio point sources often have non-trivial spectra (Herbig and 
Readhead 1992), so it is not easy to simply extrapolate from 
one frequency to another. The mean spectral index from 21 cm 
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to 1 cm is approximately a = -0.7, where S(Jy) oc v a (Cooray 
et al. 1998) but there is significant dispersion in the observed 
spectral indices (Cooray et al. 1998; Taylor et al. 2001). 

At 1 cm, there are no large scale deep surveys for point 
sources that can be used to accurately characterize the point 
source population, but there are a large number of pointed ob- 
servations toward clusters for SZ effect observations (Carlstrom 
et al. 2000) at this wavelength. The majority of observed clus- 
ters has at least one point source with a flux at 1 cm greater than 
1 mJy. When compared with the point source abundance in 
fields not containing clusters, it is clear that most point sources 
in fields with galaxy clusters must be physically associated with 
the galaxy clusters (Cooray et al. 1998), most likely the galaxy 
cluster members. 

Radio emission from galaxy cluster members has been stud- 
ied in detail (Ledlow and Owen 1996) at 21 cm {v = 1 .4 GHz), 
with the fraction of galaxies at a given radio luminosity (per 
log interval) nearly flat at radio powers (at 21 cm) below 10 248 
W Hz" 1 and falling quickly above this power. To the optical and 
radio flux limits of their survey, roughly 10% of cluster mem- 
bers showed some amount of point source emission. Assuming 
a typical spectral index a = -0.7 and our fiducial cosmology, 
the break power at 21 cm corresponds to an observed flux at 1 
cm of 7 mJy for a source at z = 0.2 and 1 mJy for a source at 
z = 0.5. This flux level is in rough agreement with the typical 
point source fluxes observed by the OVRO/BIMA SZ imaging 
experiment (Carlstrom et al. 2000; Reese et al. 2002). If 
no point source subtraction of any kind were done, these point 
sources would be a non-negligible fraction of the thermal SZ 
signal from massive clusters (~ 10 l5 h~ l Mo). 

The nearly flat probability distribution (in log flux) with point 
source flux leads to the typical galaxy cluster having the bright- 
est point source dominating the total radio flux. From an ob- 
servational standpoint this is desirable, since an SZ experiment 
only needs to remove a bright point source or two to be confi- 
dent that the measurement is not contaminated by radio emis- 
sion from galaxy cluster members. 

The expected integrated thermal SZ signal for a 10 15 Ii~ 1 Mq 
at z ~ 0.5 is on the order of 10-20 mJy (at 30 GHz) (Holder 
and Carlstrom 2001), and should scale as M 5 / 3 . This is in 
agreement with current observations (Carlstrom et al. 2000), 
but current data does not extend over a large range in mass. 
The radio point source flux should scale approximately with the 
number of galaxies, which scales as the mass (Carlberg et al. 
1996). A typical point source flux for massive galaxy clusters 
at this redshift is on the order of 1 mJy at 30 GHz and is there- 
fore roughly 5% of the total thermal SZ flux. The relative im- 
portance of point sources should therefore scale with mass as 
M" 2 / 3 , while the redshift evolution of the relative importance 
will scale as (1 +z)~ 3 for a fiat spectral index. None of these 
numbers are particularly well constrained observationally, but 
these estimates constitute our baseline model. 

4. EFFECTS OF POINT SOURCES ON THE SZ POWER 
SPECTRUM 

At some level, all clusters have radio point sources. Ideally, 
one would be able to identify and remove all point sources in 
the field of view with a beam size that is matched to the size 
of the point source and much smaller than the extent of the SZ 
signal. In such a case, the amount of SZ flux removed would be 
negligible. Not surprisingly, this is exactly the attempted strat- 
egy of choice for interferometric SZ experiments, such as the 



Ryle Telescope (Jones et al. 1993) and OVRO/BIMA (Carl- 
strom, Joy, and Grego 1996), where high angular resolution 
measurements are performed simultaneously for point source 
removal. It is not feasible to detect all point sources in a field, 
but such methods can easily remove point sources to a flux level 
below 10% of the peak SZ flux, with an effective resolution of 
roughly 10" or better. With such a strategy, the residual effect 
of point sources on the SZ power spectrum would be negligible. 

Such a strategy is not currently feasible for CMB experi- 
ments such as CBI. As primarily a CMB experiment, CBI is 
not concerned with point sources that do not contribute more 
than about 10 /iK to the rms temperature fluctuations. This 
translates into a flux threshold of a few mJy at 30 GHz. Clus- 
ters with masses near \0 u IT x Mq have a total SZ flux of a mJy 
or less (Holder and Carlstrom 2001), so a 1 mJy point source 
could be problematic. If all clusters had point sources near 1 
mJy, effectively no clusters below ~2x lO l4 h~ l Mo in figure 
1 would contribute to the anisotropy, with a somewhat reduced 
contribution from slightly larger clusters. At I = 1000 the dif- 
ference would be noticeable but at I = 4000 it might be expected 
that as much as half of the power could be missing. 

There are two leading strategies for dealing with point 
sources. One strategy is to identify possible point sources in 
catalogs at 21 cm and follow up with pointed observations at 
high angular resolution. Another strategy is to combine the data 
in a way that it is insensitive to any amount of flux coming from 
positions of known point sources, known as a constraint matrix 
approach (Bond, Jaffe, and Knox 1998). Both methods are sus- 
ceptible to point sources with "inverted" spectra, where the flux 
is higher at higher frequencies, and therefore could be missed 
in the 21 cm catalog. Such inverted sources are rare, and are 
not expected to be a dominant source of error. For extraction 
of primary CMB anisotropies both methods work quite well 
(Padin et al. 2000; Halverson et al. 2002), but for the SZ ef- 
fect, point source subtraction will remove all SZ flux within the 
point source subtraction beam area from the map at the position 
of the point source. 
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FIG. 2. — Effects of point source subtraction on SZ angular power spectrum. 
Solid line shows true power spectrum, short dashed line indicates power spec- 
tram after central 1' has been removed from all clusters. Long dashed lines 
indicate resulting power spectrum after removing central region of all clusters 
with a constraint matrix approach. This represents a worst-case scenario of all 
galaxy clusters hosting bright point sources. 
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Using the constraint matrix method, the data are combined 
such that any constant amplitude signal in the Fourier plane 
(with the phase center at the position of the point source) does 
not contribute to the measured CMB anisotropy. This will 
remove from the measured power the signal from a coinci- 
dent cluster, averaged over the data. Since the CMB measure- 
ments are sensitive primarily to scales of several arcminutes and 
higher, this will effectively remove much of the cluster signal. 
Approximately, the constraint matrix should project out of the 
data the mean cluster signal in the Fourier plane. 

As a toy model, we look at the effects of removing a point 
source from the center of all clusters. In figure 2 we show the 
effects of point source subtraction for the case of a central point 
source that has either been removed with a 1 ' (FWHM) beam or 
had power from the center of each cluster projected out of the 
data by a constraint matrix approach. In the first case, we mul- 
tiplied the SZ profiles by the appropriate beam and subtracted 
the result from the initial profile. This is correct only for an 
experiment with a reference beam that is outside the cluster, 
which is rarely strictly true, but sufficient for our purposes. For 
the second case, we calculate the mean amplitude of the clus- 
ter profile (in the Fourier domain) weighted by £ to mimic the 
effects of a uniform window function. We assume I coverage 
between £ = 500 - 4000 and assume uniform coverage of the 
Fourier plane between these values. 

The constraint matrix approach severely underestimates the 
SZ power, since it is effectively doing point source removal 
with the beam size set by the approximate angular resolution of 
the CMB measurements. Since this is larger than 1', it is to be 
expected that such an approach will remove significantly more 
power than direct subtraction with a 1 ' beam. The anisotropy 
power from galaxy clusters coincident with point sources is ef- 
fectively "nulled out" because of the relatively large synthe- 
sized beam of CBI. 
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FIG. 3. — Effect of unsubtracted correlated point sources on SZ angular 
power spectrum for three levels of point source contamination. Top to bottom: 
no contamination (solid), point source flux 2% of SZ flux for z=0.5 massive 
cluster (short dashed), point source flux 5% (long dashed), point source flux 
10% (dot-dashed). This is an approximation to the effects of no point source 
removal. Note that the power expected from point sources (which goes as £ 2 ) 
has been subtracted, leaving only the SZ power spectrum and the effects of the 
correlation. 



Direct subtraction can lead to an underestimate of the power 
of roughly 50% on scales of £ = 3000. Therefore, a detection 
of a signal of amplitude lOOfiK 2 would correspond to a true 
SZ angular power spectrum of 200fiK 2 . Improving the angu- 
lar resolution of the source subtraction helps considerably, as 
would using a reference beam within the cluster, as some of the 
SZ signal would be in both the source beam and the reference 
beam and would not be subtracted. For this reason, near future 
interferometers such as the SZA and AMI, with large spacings 
for point source removal, should not have severe problems with 
point source contamination. 

The worst-case scenario is that every point source that CBI 
projected out of the data was at the center of a galaxy clus- 
ter and that every galaxy cluster had a point source at the cen- 
ter. With nearly 100 NVSS sources per square degree and only 
about 20 clusters above 10 14 /i _1 M Q per square degree (Holder 
and Carlstrom 2001) this only requires one in five sources to be 
located in clusters. A cursory check of NVSS images of known 
Abell clusters rules out this hypothesis being correct, but it is 
not impossible that a significant fraction of galaxy clusters were 
excluded in this way. If half of the SZ signal is coming from 
clusters coincident with known point sources, CBI would be 
missing about half of the fluctuation power from the thermal 
SZ effect. 

In practice, point source subtraction will only be performed 
to a limiting flux level, meaning that many clusters will not 
have source subtraction done. In this case, the remaining point 
sources will partly fill in the SZ decrement. The extreme case 
is where none of the point sources are subtracted, leaving all 
point sources sitting in galaxy clusters, filling in the SZ decre- 
ment and reducing the fluctuation power. In fact, for the clus- 
ters contributing the bulk of the SZ signal (see Figure 1), this is 
most likely a good approximation to CBI. 

As a simple model, we assume the scalings with mass from 
section §3, with some additional constraints. We assume that all 
relevant nearby point sources are removed efficiently, and only 
include point source contamination for clusters with z > 0.2. 
This will slightly underestimate the effect of radio point sources 
for two reasons. There will invariably be some nearby clusters 
with relatively faint point sources that have escaped detection, 
and at the same time any point source subtraction, as shown 
above, can remove a significant amount of SZ flux. We also 
assume that only clusters with M > 10 l3 h~ l MQ contain radio 
sources. We assume a canonical typical point source flux of 5% 
of the total SZ flux for a 10 15 /i _1 M Q cluster at z = 0.5. 

In the Fourier domain, point source amplitudes are com- 
pletely correlated at each point, and for a central point source 
the signal will be entirely real. At any point in Fourier space, 
the combined signal from the cluster and unresolved point 
sources will be S„ e , = (T - S pt ). This will contribute to q , 
on average, S 2 et = cg- no + < S 2 pt > -if < S pl >, where 
the contribution in the absence of residual unsubtracted point 
sources. If the point sources are not concentrated at the galaxy 
cluster center (where they would be if the emission is mainly 
from the bright galaxies), but instead trace the gas we would 
expect < Spt > in the Fourier plane to have the same shape as 
f, rather than being constant. For simplicity, we assume point 
sources are strongly centrally concentrated and place them at 
the cluster centers. 

In figure 3 we show the SZ power spectrum that would be 
inferred. The power spectrum due to the point sources alone 
has been subtracted. Unsubtracted point sources cause a signif- 
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icant underestimate of the SZ angular power spectrum. Assum- 
ing a larger contamination of 10% of the SZ flux for a z = 0.5 
10 15 /z _1 M Q eliminates nearly all of the SZ power, while drop- 
ping the contamination to only 2% reduces the missing power 
to only about 20%. 

5. DISCUSSION 

Point source contamination makes interpretation of detection 
of fluctuations in the CMB from the thermal SZ effect very dif- 
ficult. To compare these fluctuations to predictions from ei- 
ther semi-analytic modeling or numerical simulations, some 
treatment of the effects of point source subtraction (or non- 
subtraction) is required. This would require a recipe for the 
cluster radio galaxy populations. Current predictions of the 
thermal SZ power spectrum are almost certainly overestimat- 
ing the thermal SZ power that CBI could observe. 

The measured power at 30 GHz at high multipoles could be 
less than 30% of the true SZ power, but more likely is measur- 
ing roughly 50-75%. Specifically, the recent tentative report of 
temperature fluctuations at high multipoles of close to 500/j.K 2 
(Mason et al. 2002), if due to the thermal SZ effect would 
indicate a true signal on the sky of roug hly 1000^ 2 . A tem- 
perature rms of 25 [iK represents a bit of a challenge for theo- 
retical models (Bond et al. 2002; Komatsu and Seljak 2002), 
so a true signal greater than 30^ would require a major re- 
thinking of the physics of galaxy clusters and/or some fine tun- 
ing of cosmological parameters. The simplest way to increase 
the expected SZ power is to increase a%. A true SZ signal of 
more than 30/iK would suggest as > 1.1 (Komatsu and Seljak 
2002), a value that is not preferred by current CMB data (Bond 
et al. 2002) but is not ruled out. From figure 3 of Bond et al. 
(2002) such a high value could most easily be accommodated 
if the Hubble constant were significantly lower than the value 
suggested by HST measurements (Freedman et al. 2001). 

Alternatively, if the measured power at high £ really is a mea- 
surement of the thermal SZ effect, this would be evidence that 
low-mass clusters at z ~ 0.5 are remarkably devoid of bright 
radio sources. This suggests that upcoming SZ surveys at 30 
GHz (SZA) or 15 GHz (AMI) should have surprisingly clear 
extragalactic skies. 

The statistics of point sources at high radio frequencies are 
very poorly constrained, making detailed predictions of point 



source contamination difficult. We have adopted an approach 
to modeling point source contamination of SZ signal that is em- 
pirically motivated, with only the relative importance of point 
sources to SZ signal at a single mass scale as a free parameter. 
Upcoming SZ/CMB experiments with high angular resolution, 
specifically for the purpose of point source detection, such as 
AMI and SZA will provide a wealth of information on radio 
point sources at these frequencies, while at the same time pro- 
viding valuable information on the SZ effect from galaxy clus- 
ters. 

Single frequency measurements of the fluctuations due to the 
thermal SZ effect at frequencies below ^90 GHz that do not 
subtract point sources with a small beam to a fairly low flux 
level will be contaminated at a largely unknown but almost cer- 
tainly significant level. Higher frequency measurements could 
have similar problems from dusty starburst galaxies, but most of 
these sources are expected to not be associated with the galaxy 
cluster members. While lensing effects lead to an enhancement 
of the confusion noise (Blain 1998), it does not lead, on aver- 
age to an increased average flux. 

Experiments with multiple frequencies will be required for 
a robust determination of the amplitude of the thermal SZ sig- 
nal, and the strong correlations between radio point sources and 
galaxy clusters and submm point sources and galaxy clusters 
(primarily due to gravitational lensing) will require careful at- 
tention. 

The general problem of correlation between secondary 
anisotropics and foregrounds will be increasingly important. 
The thermal SZ effect may be the most significant example, 
but correlations between the lensing of the CMB, for example, 
and radio and submm point sources will reduce the expected 
signal and/or modify the noise properties of any attempted re- 
constructions. Clearly, a better understanding of the covariance 
between various foregrounds and backgrounds and secondary 
anisotropics of interest will be required if such signals are to 
be used as useful tests of our understanding of cosmology and 
structure formation. 

This work was supported by a W.M. Keck Fellowship at the 
IAS. I would like to thank Erik Reese for many useful conver- 
sations. 
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